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Abstract. We explain the peak spacing in the power 
spectra of millisecond variability detected by the Rossi X- 
ray Timing Explorer (RXTE) in the X-ray emission from 
the LMXBs (Sco X-1, 4U1728-34, 4U1608-522, 4U1636- 
536, etc) in terms of the rotational splitting of an intrinsic 
frequency (which is of order of the local Keplerian fre- 
quency) caused by an accretion disk. We calculate this 
effect and demonstrate that there is a striking agreement 
with the observational data. We show that the observed 
discrete frequencies ranging from 200 to 1200 Hz can be 
described by a whole set of overtones. For higher overtones 
(lower frequencies, ^ 100 Hz) the ratio between frequen- 
cies is determined by the quantum numbers alone. We sug- 
gest that a similar phenomena should be observed in Black 
Hole (BH) systems for which the QPO (quasi-periodic os- 
cillation) frequency should be inversely proportional to 
the mass of the compact object. For BH systems the char- 
acteristic frequency of oscillations should therefore be a 
factor of 5-10 less than for a neutron star system. 

Key words: Accretion: disk-binaries: LMXBs -X-Rays: 
QPOs 



1. Introduction 

The launch of RXTE opened up a new era in the study of 
quasi-periodic oscillations (QPOs). Recently, kHz QPOs 
have been discovered in the persistent fluxes of eight 
LMXBs: Sco X-1 (van der Klis et al. 1996a), 4U1728-43 
(Strohmayer et al. 1996a), 4U 1608-52 (Berger et al. 1996), 
4U 1636-53 (Zhang et al. 1996), 4U 0614-1-091 (Ford et 
al. 1996), 4U 1735-44 (Wijnands et al. 1996), 4U 1820-30 
(Smale et al. 1996), and GX5-1 (van der Klis et al. 1996). 
In addition, episodic and nearly coherent oscillations have 
been discovered during several type-I X-ray bursts from 
KS 1731-260 with a frequency of 363 Hz (Strohmayer et 
al. 1996a), during one type-I burst from KS 1731-260 with 



a frequency of 524 Hz (Morgan & Smith 1996). In ad- 
dition during 3 bursts from vicinity of GRO J 1744-28, 
oscillations with a frequency of 589 Hz (Strohmayer et 
al. 1996a), and during 4 type-I bursts from 4U 1636-53 
with a frequency of 581 Hz (Zhang et al. 1996) were seen. 
The presence of a pair of kHz QPOs is characteristic of 
almost all these observations. The centroid frequency of 
these QPOs ranges from 400 (4U 0614-1-091, Ford et al. 
1996) to 1171 Hz (4U 1636-53, van der Klis et al. 1996c). 
Also, it is worthwhile to point out that van der Klis et al. 
(1997) reported that they detected, simultaneously with 
the kHz QPO, two additional peaks near 40 and 90 Hz. 

In this Letter we shall demonstrate that the effect of 
a rotational splitting (caused by the Keplerian accretion 
disk) of oscillation frequencies that has been ignored in 
all previous theoretical studies of QPOs may be important 
for our understanding of a phenomenon of QPOs in binary 
systems with neutron stars and black holes. Our conclu- 
sion about relevance of this effect to the oscillations in the 
innermost region of accretion disk around a compact ob- 
ject is robust and does not depend on a particular model 
of a disk. If the QPO phenomenon is associated with such 
oscillations, then we must observe the rotationally split 
frequencies. The detailed study of this effect may provide 
us with very important information about the physics of 
accretion disks. 



2. The Effect of Rotational SpHtting by an Accre- 
tion Disk 
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It is well known (see e.g. Unno et al. 1979, and references 
therein) that in a rotating star the non-radial oscillations 
are split: (Tkim = o'ki + m51*Cki, where dki is the frequency 
in the nonrotating star, tJkim is the frequency in a rotating 
star (in the corotating frame), m is an azimuthal number, 
fl^ is the stellar rotation frequency, and C^i is expressed in 
terms of an integral that depends on the stellar structure 
and on the eigenfunction. Thus, in a rotating star, a non- 
radial oscillation drifts at a rate mJl^tCki relative to a fixed 
longitude of the star. The effect of rotational splitting is 
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due to the Coriolis force and is similar to the Zecman effect 
in a magnetic field. 

The oscillations of a rotating disk must split as well. 
The rotational splitting of the disk oscillations is partic- 
ularly important for those modes whose frequencies are 
comparable to the local rotation frequency of a disk rela- 
tive to a distant observer. The detailed theory for a star 
is given by Ledoux & Walraven (1958) and by Unno et al. 
(1979). Here we present calculations of the splitting effect 
for an oscillation frequency due to the Coriolis force in a 
disk geometry. 

The general formula for the correction to the oscil- 
lation frequency of a star/disk produced by the Coriolis 
force reads (cf. Unno et al. 1979, equation [18.28]) 



(0) >.(0)* 



Jv^k,m^l 



k,m 



dV 



(1) 



Here n is the local angular velocity in a disk in the lab- 
oratory frame of reference, ^^^^s^{t,T, ip,7,) is the (k,m)- 
displacement component, and the integration is over the 
volume of a disk/ring-like configuration. 

Let us consider the oscillations localized in some annu- 
lus of a disk at radius R. We assume that the amplitude 
of these oscillations decreases toward the innermost edge 
of a disk. Then, the component of a displacement vector 
can be calculated by using the (k, m)-harmonics of a com- 
plete set of eigenfunctions for a disk/ring-like configura- 
tion {Uk,m = ruk,m(r, z, t)}, where r, tp, z are cylindrical 
coordinates, and t is the time in a laboratory frame of 
reference. In equation (1) the integration has to be taken 
over the volume of a ring of radius R, width AR, and 
half-thickness H. The vectors ^|||\jj-)(t, r, (p, z) are usually 
presented in such a form where their components are pro- 
portional to the corresponding components of a gradient of 
a potential Uk,m — ruk In cylindrical coordinates they 
can be written as (cf. Unno et al. 1979, equation [18.28]) 
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The functions Uk.m satisfy the free-boundary conditions at 
z = ±H (see e.g. Morse & Feshbach 1953): 
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Also, these functions must be symmetric about the disk 
plane 



Uk, 



i(t,Z,</?) = Ukan(t, -Z, (^), 



where < z < H. 

In the case of uniform rotation, the temporal and az- 
imuthal dependences of eigenfunctions Uk,m are given by 
exp[i(m(p — riot)], where fio is the angular eigenfrequency 



of the oscillation (see Unno et al. 1979), thus, the compo- 
nent Uk m can be written as 



Uk,m = e'™''' cos(7rkz/H)e 



-iOnt 



(5) 



where m — ±1, ±2, and k — 1,2, ... . 

Note that the functions Uk,m are eigenfunctions of a 
two-dimensional Laplace operator. They satisfy the free- 
boundary and symmetry conditions over a z-coordinate 
(see equations [3], [4]) and the periodicity condition over 
a (^-coordinate. It is important that these eigenfunctions 
form a complete set of functions (a fundamental system) 
which can be used as a basis for the expansion (see e.g. 
Morse & Feshbach 1953) of any arbitrary displacement 

The numerator and denominator in formula (1) are 
given by 



Ii = -27rHRAR(20m^rCv): 
and 

h = 2^HRAR[g + {i.R/iifTT\^ + m'^^l]. 



(6) 



(7) 



respectively. Here ^r, £.<p, and are the average- weighted 
displacements for a disk/ring-like configuration in radial, 
azimuthal, and vertical directions, respectively. To illus- 
trate the effect of rotational splitting for a disk/ring-like 
configuration we may justifiably assume that ~ ^<p, then 
we arrive at 



.(1) 



2mn 



(8) 



where s(R/H) — (CzR/CiH)^ is a function determined by 
the vertical structure of a disk. For more or less realistic 
structure of a disk s 1 and depends on the ratio R/H, 
which is either ^ 1 (near a compact object) or ^ 1 (far 
away from a compact object). We must note that the re- 
lationship between ^r, and depends also on the az- 
imuthal and vertical mode numbers. For example, consider 
an individual oscillation mode with displacement Cfk^m) 



»(k,m) 

an incompressible fluid. Then V • ^'"^^-j = 0, and we get a 
linear partial differential equation that determines a rela- 
tionship between ^r, and ^z- This relationship, accord- 
ing to expression (2), should depend on the azimuthal and 
vertical mode numbers, m and k, respectively. 

Thus, the frequency of the oscillations dete cted by a 
distant observer is given by (cf. Unno et al. 1979 , equation 
[18.33]) 



(4) fik,m = f^o + mO + o-(i) = 



fio + mO 1 



STT^k^ + m^ + 1 / ' 



(9) 



where m and k are the azimuthal and vertical mode num- 
bers, respectively. 
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3. Explanation of QPO Frequencies in Terms of 
Rotational Splitting 

To estimate the frequency splits, we may adopt CIq = Q, 
which is a very good approximation for a nearly Keplerian 
accretion disk. In this case, for the lowest-order modes 
with m = 0, —1 and k = 1, 2, 3, 4, and 5 the oscillation 
frequencies can be calculated as 

Oi _i = 217o/(s7r2 + 2), 
02,-1 = 2f}o/(s47r2 + 2), 
^3,-1 = 2r!o/(s9^2 + 2), 
^4,-1 = 2no/(sl67r2 + 2), 

and 

fis-i = 2f]o/(s257r2 + 2), 

respectively. 

For m = — 2 and k = 1 we have 

|fil,_2| = rio[l- 4/(877^ + 5)]. 

Let us assume that uq = Q,q/2-k = 1200 Hz. Then, for 

the most plausible range of s w 0.7—1 we get the following 
frequencies (arranged in ascending order, in Hz): 



and 
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These frequencies perfectly match the observed QPO fre- 
quencies in LMXBs. It is important that for the modes 
with a low azimuthal number m = — 1 and high k- numbers 
(high overtones in the vertical direction) the ratio be- 
tween frequencies does not depend on the function s, 
rather it is solely determined by a quantum number k: 
i/k,_i/!/k+i,_i w (k-l- l)^/k^ (k > 2). The frequencies of 
modes with m = —2 and k > 2 concentrate to 1^1^.-2 ~ vq. 
It must be pointed out that the modes with m = 1, 2, 3, ... 
and k > 1 have frequencies concentrating to the frequency 
(m + 1)vq. Perhaps, these modes are not allowed or, for 
some reason, are not excited. In this regard, it would be 
interesting to investigate the selection rule for the oscilla- 
tions in the innermost part of accretion disk. 

It is important to note that a contribution of each in- 
dividual (k, m)— component to the power spectrum of os- 
cillations is determined by the smoothness of a function 
characterizing the perturbed surface. It is well-known from 
the Fourier analysis that the power of the /-component is 
~ , where N is the order of highest existing derivative 



of a function describing the perturbed surface with respect 
to the corresponding coordinate. In the case of the oscil- 
lations of the innermost part of accretion disk, / = m or 
Z = k for the azimuthal or vertical perturbations respec- 
tively. It is very likely that the azimuthal perturbations 
are smooth, so that we may expect the presence of only 
the lowest order modes with |m| = 0,1,2 in the power 
spectrum. On the contrary, the vertical perturbations are 
essentially discontinuous, and their Fourier spectra should 
contain a large number of modes with k > 1, with the con- 
tribution of each of these modes to the power spectrum 

~k-2. 

The occurrence of the localized oscillations in the disk 

may indicate that there is a local release of (e.g. gravita- 
tional) energy in addition to the energy of viscous dissi- 
pation. If so, then the local plasma temperature should 
be T ^ 10 — 20 keV, which is much higher than the 
temperature in a standard viscous disk (see Shakura & 
Sunyaev 1973). The possibility of such oscillations, their 
energetics, high quality Q, and substantial modulation of 
the observed X-ray luminosity due to these oscillations, 
are discussed by Titarchuk et al. (1997) in the context of 
the effect of a centrifugal barrier which is thought to take 
place in a boundary region surrounding a neutron star. 

4. Summary 

We suggest that in the QPO phenomenon wc observe the 
effect of a rotational splitting of the oscillation frequency. 
For a Keplerian disk this effect should be most pronounced 
because the characteristic frequency of the oscillations of a 
centrigugal-barrier region is of the order of the Keplerian 
frequency. 

We conclude that the observed three frequencies rang- 
ing from 300 to 1170 Hz can be naturally interpreted in 
terms of a rotational splitting of the main frequency: the 
lowest frequency, ^ 200-300 Hz, corresponds to a mode 
with m = — 1 and k = 1; the higher frequency, ~ 800 Hz, 
corresponds to a mode with m = — 2 and k = 1, and the 
highest frequency is the main frequency of the oscillations. 
The issue of the excitation of these modes and emission 
mechanism will be addressed elsewhere. 
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